Probing photon energy statistics of an emitter using photon correlations 



A. Trichet 1 *, S. Bounouar 1 ' 2 * 
1 CEA-CNRS-UJF group 'Nanophysique et Semiconducteurs', 
Institut Neel, CNRS - Universite Joseph Fourier, 38042 Grenoble, France, 
2 CEA-CNRS-UJF group 'Nanophysique et Semiconducteurs', CEA/INAC/SP2M, 38054 Grenoble, France 

We investigate the possibility of measuring the homogeneous and inhomogeneous contribution 
to the linewidth of a spectrally diffusing single photon emitter using a simple photon correlation 
spectroscopy method (PCS). The photon energy statistics of the homogeneous line (poissonian) and 
of the spectral diffusion (first order markovian) being of different natures, they act differently on 
the half- line autocorrelation function (HLAF). We model here their effects and show it is possible to 
extricate them, opening the opportunity to determine separately the homogeneous linewidth, and 
the spectral diffusion amplitude. 
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I. INTRODUCTION 



It is known since the pioneer works of Handbury-Brown 
and Twiss on the correlation of an optical field [1], and of 
Kimble, Dangennais and Mandel on the correlation of a 
quantum state of the light |2J , that it is possible to probe 
the emission statistics of an emitter by photon correla- 
tions. In particular Kimble et. al showed that for a two 
levels system, the autocorrelation function of the emit- 
ted optical field increases with the delay r |3J. This is a 
characteristic of a single photon emission. In condensed 
matter systems, the energies of the emitted single pho- 
tons can fluctuate. The interaction of the emitter with 
the numerous degrees of freedom of the surrounding en- 
vironment leads to random fluctuations of the emission 
energy. Since quantum computation requires indistin- 
guishable photons [4] in order to proceed to efficient two- 
photons interferences [5], this is a severe limitation in 
the use of solid state systems for such operations. This 
dephasing phenomen is the result of several different pro- 
cesses. From the energy statistics point of view, one can 
divide them in two categories: Poissonian processes, such 
as exciton-phonon coupling, randomly affect the energy 
of the emitted photons. They are "memory less" pro- 
cesses as the energy of the emitted photon at time t does 
not depend on the energy of the emitted photon at time 
t — r whatever the value of r. We define the energy distri- 
bution of these uncorrelated random processes as the ho- 
mogeneous linewidth of the emitter. On the other hand, 
a correlated process can shift the emitter energy. Fluc- 
tuations of the electronic environment are also a source 
of dephasing. This random Stark-shift of the emission 
energy is called spectral diffusion [6]. In opposition with 
phonon coupling, spectral diffusion is time-correlated and 
it is is a reasonable approximation to consider it as a first 
order Markovian process with a correlation time r c [ZJ|8]. 
Unlike the poissonian process, which is independent of its 
history, the first order markovian process is influenced by 
its immediate or most recent past. Thus, the energy of 
the emitted photon shifted by the spectral diffusion de- 
pends on the energy of the previously emitted photon. In 



a previous publication [9j, it was shown that, by taking 
the Handbury-Brown and Twiss setup and by adding de- 
tection energy conditions, such as selecting only one half 
of the emission line, one has access to the subnanosecond 
correlation time of the spectral diffusion r c |9J. In this 
letter, we show that this measurement technique is even 
more powerful since the resulting half line autocorrela- 
tion function (HLAF) brings informations on the rela- 
tive contribution of the correlated and uncorrelated pro- 
cesses. Thus, such a slight modification of the original 
HBT setup, which was providing emission statistics of 
the emitter, gives as well access to the energy statistics 
of the emitter. In a first part, we describe the main prin- 
ciples of the experiment. Then, after showing that it is 
always possible to separate the emission statistics and 
the energy statistics contributions on the HLAF, we de- 
rive analytically the HLAF of a spectrally diffusing single 
photon emitter with a finite homogeneous linewidth. Af- 
ter discussing its properties, we confirm the validity of 
our model with a Monte Carlo simulation, building nu- 
merically the HLAF of a spectrally diffusing poissonian 
emitter and simulating the influence of a finite homo- 
geneous linewidth. To finish, we highlight the possible 
direct application of this theoretical result: the sepa- 
rate determinations of the homogeneous linewidth and 
the spectral diffusion amplitude. 



II. PRINCIPLE OF THE MEASUREMENT 

The principle of the PCS experiment is to convert fluc- 
tuations of the emitter energy in intensity fluctuations. 
In fig. [I] a) we consider the situation where all the pho- 
tons are collected and sent to the HBT setup whatever 
their energies (photons comming from the right side of 
the line are coloured in red, from the left side in blue). 
The scheme in fig. [I] b) represents the single photon 
emission with photons arriving one by one with a char- 
acteristic time defined by the emitter radiative lifetime. 
The HBT setup allows to measure the intensity first order 
correlation function |10| : 
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Figure 1. a) spectrum of a spectrally diffusing emitter, b) 
photons statistic of the emitter in the HBT experiment c) 
Autocorrelation function of the single photon emitter convo- 
luted with the temporal resolution. 

This function is plotted in fig. [l]c). The dip at zero 
delay is the signature of single photon emission. It is 
non-null at zero delay only because of the convolution 
with time resolution of the photodetectors. 

We introduce now a detection energy condition such as 
we only detect photons emitted from the high energy side 
of the spectrum called right side (R side) in fig. [2] There- 
fore, we can discriminate the photons emitted from the 
high energy side of the spectrum (ie. R side, red spikes 
in fig. [2] b) ) from the photons emitted in the low energy 
side (L side, blue spikes in fig. [2]b) ). In the fig. [2] b) , 
bunches of single photons (red spikes) emitted from the R 
side are observed because the spectral diffusion is a time- 
correlated process (ie. with memory). On the HLAF (fig. 
[5]c)), it results in a bunching combined with a zero de- 
lay dip due to the subpoissonian nature of the emission. 
If we consider that the emitter homogeneous linewidth is 
infinitely small, the HLAF can be easily calculated with a 
rate equation model by considering a split two-levels sys- 
tem with a probability for the emitter energy to switch 
from one side to the other one. 



HLAF(r) = g 2 (R,R,r) = [l+exp(-j c r)} [l-exp(-(r+j)r 

(2) 

with 7 C = the correlation rate, r the pump rate and 
7 the spontaneous emission rate. 

It is the product of two expressions. Since [1 — 
exp( — (r + 7)t)] describes the single photon behavior of 
the emitter, [1 + exp(— 7 c t)] describes the spectral diffu- 
sion phenomenon and its influence on the ability of the 
emitter to send photons in the right spectral window. 
The spectral diffusion term is equal to 2 for the null de- 
lay. 
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Figure 2. a) spectrum of a spectrally diffusing emitter with an 
infinitely small linewidth: only the right spectral window is 
selected for detection, b) photons statistics of the emitter in 
the PCS experiment: only red photons are detected c) HLAF 
in the case of an infinitely narrow linewidth 



We now consider a finite homogeneous linewidth (fig. 
|3| . The mean value of the homogeneous lineshape tran- 
sits from one side of the spectral window to the other 
with the time constant r c . However, even if it is centered 
in the left side, a photon can be emitted in the right side 
and vice versa. The consequences on the photons statis- 
tics are represented in fig. [3]b). To the bunching statis- 
tics due to the time correlated spectral diffusion process 
it has to be added an uncorrelated random distribution 
of the photons on both spectral windows, degrading the 
photons "bunching". The consequence on the HLAF is 
a diminution of the bunching contrast, as the energies of 
photons are less time correlated (fig. [3]c)). The impor- 
tance of this contrast is linked to the relative importance 
of the homogeneous linewidth compared to the spectral 
diffusion amplitude. To understand this dependence, we 
calculate in the following the HLAF corresponding to a 
spectrally diffusing single photon emitter with a finite 
homogeneous linewidth. 



III. ENERGY AND EMISSION SEPARATION 

The aim of this section is the obtention of a preliminary 
] result, which will simplify the calculations leading to the 
derivation of the HLAF for a non-negligible homogeneous 
linewidth. 

We propose to show that the HLAF is the product of 
the spectral diffusion part (carrying information on the 
correlation of the emitter energy) and the single photon 
part (carrying information on the correlation of the emis- 
sion). This result only holds when emission and energy 
position of the line are independent. The emitter is seen 
as a two-level system (TLS) of fluctuating energy. 

Let us define the following ensembles: 
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Figure 3. a) spectrum of a spectrally diffusing emitter with 
a finite linewidth: only the right spectral window is selected 
for detection, b) photons statistics of the emitter in the PCS 
experiment: only red photons are detected c) HLAF in the 
case of a non-negligible homogeneous linewidth 



-7 £ : ensemble of events such as a photon is detected 
at time t. 

- 7g m : ensemble of events such as a photon is emitted 
by the TLS. 

- r )\ n \ ensemble of events such as the TLS energy is in 
the right spectral window. 

The emission of a photon by the two level system (TLS) 
and its energy situated in the right spectral window are 
two independent events, so: 



(3) 



We go back to the expression of autocorrelation func- 
tion in term of probabilities: 



,2,x _ P(Y I 7°) 



p(l T ) 



(4) 



To detect a photon at time t, one needs to have a 
photon emitted by the TLS at time t and to have the 
TLS energy in the right spectral window at time t: 



lem ^ 7in 



Therefore, 



p(Y | 7°) = Phlm Yin I lem H lL) 

Using ([3| in we have: 



(5) 
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We can now find the general expression, separating 
the TLS emission and energy correlation functions in the 
>HLAF. 
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This result will make our future calculations of HLAF 
easier, since the emission part ^ m (r) will always be the 
same and we will only have to focus on the derivation of 
9in( T )i tne second order correlation function of the TLS 
energy position. 



IV. HALF-LINE AUTOCORRELATION 
FUNCTION 

Since the mean TLS energy fluctuates from the right 
side to the other with a correlated statistics, but can also 
be randomly, and with no time-correlation, distributed 
in energy along the lorentzian profile imposed by the ho- 
mogeneous linewidth, The total energy £^of the TLS at 
time t is the sum of two random variables: 



P 



(9) 



where // is the energy position of the center of the 
homogeneous lineshape at time t. It is distributed along 
the gaussian distribution of the fluctuations jT2j[13] . This 
random variable is time correlated. It means that the 
value of the random variable at time t + r is influenced 
by its value at time t and : 



< M y +T > <//><M t+T > 



(10) 



e l is the energy shift due to the homogeneous linewidth 
at time t. We assume it has a lorentzian distribution 
centered on /A This variable is described by a poissonian 
process and is not time-correlated. 



{£ t £ t +1 = {£ t ){£ t +1 = {£ o ){£l = {£ ry. 



(11) 



The spectral window of detection is defined by the en- 
ergy interval I n . In this letter, we consider, I n = [0, +oo[, 
ie. corresponding to the R side. 

The energy of the TLS is in the right spectral window 
at time t when E t el n ^ so when: 



Therefore, we have: 



■{(et + SWn)} 
{(s° + /i°)e(/ n )} 



(12) 



P ( 7 (r) | 7 (0)) = P (Yem I 7em H lL)p{lln I Tern H 7"J 

= Pinlm I 7em)p(7* T n I lin) (7) 



We also define the following ensembles : 
-/4 n : ensemble of \i such as /i £ e/ n . 
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-Mont : ensemble of /i such as // ^ I n 
We found out in the previous subsection that we can 
always write down the autocorrelation function as: 



9\t) = 

with 9 Ur) = 



2 

9 em 



Pilin I 7?„) = p{{{e T + p T )e(I n )} | 7 °J (13) 

We use here the law of total probability which asserts 
that for an ensemble A and its partition {Ai} such as 
A = J2 i A i , 



P(B) = J2P(B\A i )P(A i ) 

i 

So, eq. ( fi~3| ) becomes: 

Pilln I ll) = P({(£ T + P T Wn)} | ML I iDpipln I 7°„) 
+ P ({(^ + » T Wn)} I Mouf I ll)piPlut I 7°„) 

(14) 

({(e r + fi T )e(I n )} | /i[ n ), the ensemble of values taken 
by e and \i such as [e + /i)e(J n ) at time r, knowing that 
/ie/ n at time r is independent from the possible values 
taken by e at time 0. Indeed, e is a poissonnian random 
process. Moreover, the condition on the \i value is already 
fixed at time r in ({(s T + /i r )e(/ n )} | Mln)> a condition 
on its value at time does not change the ensemble. 

So we can conclude that the ensemble 
({(e T + fi T )e(I n )} | fjLj n ) and j® n are independent 
ensembles, thus: 

P ({(e T + fl T Wn)} I ML I 7°>(mL I 7°J 

= p({(e T + M T )e(/n)} | plnMPln I 7?J (15) 

and 

Pilln I 7to) = P({(^ +p T Wn)} I /4>(/4> I 7°„) 

+ P ({(e T + fWn)} I PlutlPiPlut I 7°„) (16) 

For any ensembles A and B, one have the basic rela- 
tionship: 



P(A | B) = P(B | A) 



P(A) 
P(B) 



so (16) becomes: 



Pilfin I Pin) = P((7°n I Mi J I Pin)p(Pin I ML) 
+ P((7°n I Mont) I ViJpilfiut I ML) 
= P(7*°n I Pin)p(Pin I Mln) 
+ P(7z°n I MoutMMout I Pin) 

= p(iL\pI)p(pL\pI)^^ 

p(Pjn) 

^p(iL\pLt)p(pL\pLt)^^ 



(18) 



The same way, 



Pill I Pint) = Pill I I pI)4tPs 

PKPout) 

+ Pill I pLMplut I mD^H (is) 

P\Pout) 

We inject the two last equations in ( fl7| ): 

Pilin I 7in) = Pilin I Pin)p(lin I Pin)p(Pin I Min) j 
+ P(7,n I Mln)p(7?n I /oJ^^n I Mo"nt) 
+ P(7[n I MLt)P(7?n I tflnWtimt I 
+ P(7ln I MSuXtSi I PluMPlut I 

(20) 

P(Pin)> p(Pout)& nd p{lin) are respectively the proba- 
bility to have /i G I n , /i ^ I n: and {(£ r + /i r )e(/ n )} at 
null time. Thus: 



P(Pin) = p(Pout) = Pi,! in) 



1 



and (20) becomes 



Pilln I Jin) = Pilln I Pin)p(lin I Pin)p(Pin I Min) 
+ P(7z T n I Pin)phin I mLMmL I Mont) 
+ P(7in I Mont)p(7z°n I Pin) ] P^Plut I Min) 
+ p{lin I P T out)p(lin I Pou^PiPlut I Mont) 



(21) 



M(7ln I ll) = P({(e T + P T W n )} | /ilJ/^n I Mln)^y 

+ p(i(eT + lS)<In)} I P T out)p(lin I ^ont)^y (17) 

We apply again the law of total probability on the 
terms p(^ in \ 

Pin ) and p(7° n | 

Pout)- 



The probability for the TLS to be in a the right spec- 
tral window, knowing that it was in it at t=0, is the sum 
of four probabilities, describing the four possible config- 
urations for the occurrence of a coincidence: 

-For the first term, the homogeneous line is in the R 
spectral window at t=0 and also at t=r. 

-For the second, the homogeneous line is out of the R 
spectral window at t=0 and in at t=r 



-For the third, the homogeneous line is in the R spec- 
tral window at t=0 and out at t=r. 

-For the fourth term, the homogeneous line is out of 
the R spectral window at t=0 and also at t=r. 

Let us evaluate p(jj n | Pj n ( out ))i the probability that 
the TLS energy is in the right spectral window at time 
r when the center of the homogeneous line is in (out 
of) the right spectral window. Despite the fact that the 
homogeneous linewidth is of relative complex shape 
we approximate it as a lorentzian centered on E — \i. In 
accordance with the Kubo-Anderson model, the energy 
distribution of the fluctuations is gaussian [12j 03] • We 
have: 



Pin(out)) 



P(Pin(o 



it)> 



2 p(lin n Pin(out)) 
(22) 



We consider a Lorentzian centered in the right spec- 
tral window (ie. \ieli n \ and calculate the probability to 
measure a photon emitted in the right spectral window. 



R side Gaussian distribution 
of the line energies 
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Figure 4. Lorentzian centered on energy \i and gaussian distri- 
bution of the possible \i positions. This distribution is cut by 
the spectral selection of the right part of the inhomogeneous 
line. It defines the area of the lorentzian where photons can 
be detected. 

For a given lorentzian centered in /i, the probability for 
the TLS to emit a photon in the right spectral window is 
the ratio between the lorentzian area in the right spectral 
window and its total area: 



with 



f°°lor(<T,E- n)dE 
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JEeI iri 



lor{a, E 



H)dE 
(23) 



We have then to consider the probability for the 
lorentzian to be centered at the energy E between \i and 
ji + d\i. It is obtained in multiplying the previous ex- 
pression by the probability density imposed by the gaus- 
sian distribution Gauss(Y>, fijdfi and integrating over all 
fx such as fieli n . 



P{lin ^ Pin(out)) — 
Lzlin(out) GaUSS ( E ' p) fEeI in ( l0r ( a ' E ~ P) dE ) d P 

/ °° Gauss(Yi, n) J °° (lor(a J E — ji)dE)dfi 



(24) 



with Gauss^, fi) 



= exp 



(2^)' 



E a V27r~~ jr V 22J « J ' ~* 2V2ln2 

the standard deviation and E the full width at half max- 
imum of the gaussian function. 

This is the probability for the TLS to be in the right 
spectral window and the lorentzian centered in the right 
(left) spectral window at any time, so it does not depend 
on the delay r and are only functions of the homogeneous 
linewidth a and the gaussian standard deviation E5: 

pin T in^P T in{out)) = Phin^Pin(out)) = OL in{out) ((J,Y>) (25) 

We rewrite the eq. [2l]with these coefficients: 

p(lin I Jin) = ^[ a inP(Pin I Pin) 

+ Qtin a outp(Pin I Pout) 
+ aout^inp(p T out I Pin) 
tPiPlut I Pout)] 



(26) 



This probability depends on one hand on the coef- 
ficients Oii n and a out which depend themselves on the 
linewidth of the homogeneous lorentzian line and of the 
gaussian distribution. They take into account the ran- 
dom uncorrelated statistics introduced by the finite ho- 
mogeneous linewidth. On the other hand, eq. ( [26] ) 
depends on the probabilities p(fij n | Mi n )> p(Pin I 
P Q out)APlut I Pin), PiPlut I Pout) which are describ- 
ing how the mean value of the homogeneous linewidth 
is transiting from one spectral window to the other and 
take into account the time-correlated part of the statis- 
tics. All these probabilities can be determined analyti- 
cally with the model of the infinitely sharp homogeneous 
linewidth presented in the reference |9J , and its associated 
rate equations. 

Because of the configuration chosen in this letter (ie. 
left side versus right side of the spectrum), the proba- 
bility 7# for the mean value of the homogeneous line to 
jump from the right spectral window to the left is the 
same as the probability for the homogeneous line to 
make the opposite move, thus: 



6 



Then, we can derive: 

p(Pin I Vin) = — +(1-— )exp(-j c r) = ]- (1 + exp(-j c r)) 
7c 7c ^ 

(28) 



p(Pin I ML) = yl 1 - exp(-j c r)} 
= i[l - exp(-7 c r)] 



Pitfout I ML) = ol 1 + e ^(-7c 



We put all these expressions in eq. ( 26 ) : 



(29) 



PiYin I 7z°n) = ^ ( a in + °W) + 

+ (^(«in + a lut) ~ a in a out ).exp(-j c r)} (30) 

To obtain the correlation function, we divide this ex- 
pression by p( 7 ^ | 7°J = 4[|(a? n + o£ wt ) + a in <w], 
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Figure 5. Bunching factor /3 as a function of the fluctuation 
amplitude E, and the homogeneous linewidth a 



so, P —> 1, which is its expected value in the case of 
the infinitely sharp linewidth. 

When the homogeneous linewidth becomes larger than 
the standard deviation of the fluctuation, a ^> E: 



— » 



4 W = 1 + 



( 2 ( a in + a ont) Qt-in&out) 



exp(-<y c r) (31) 



finally, rearranging terms and multiplying the emission 
part <7em( r )j one nas: 



£ 2 (i?,i?,T) = [l + 



-.exp(-7 c r) 



* [1 - exp(-(r + 7)r)] (32) 
Thus, the HLAF is on the form: 
g 2 (R, R, t) = [l+/3.exp(- 7c r)][l-exp(-(r+7)r)] (33) 
with, 



/3 



(a* 



(a in + ow) 2 



(34) 



We can make two remarks about this factor: 

- The bunching factor f3 is exclusively ruled by the 
uncorrelated random statistics of the finite homogeneous 
linewidth (expressed by coefficients c^ n and a out ), and is 
a function of a and E. As shown in fig. [5j the bunching 
is important (red color) for small a and large fluctuation 
amplitudes E. 

- The behavior of j3 at the limits is interesting: when 
the homogeneous linewidth tends to 0, ie. a <C E, we 
have: 



2 




and /3 —> 0, the bunching is collapsing. This situation 
corresponds to a case in which the poissonian statistics 
of the homogeneous linewidth takes over the correlated 
statistics of the spectral diffusion. 

The HLAF bears simultaneously the signature of the 
subpoissonian emission statistics (zero delay dip), the 
correlated spectral diffusion energy statistics (bunching) , 
and the poissonian energy statistics of the homogeneous 
linewidth (degradation of the bunching), the two last sig- 
natures having opposing effects. 



V. MONTE-CARLO SIMULATION 



In order to confirm the analytical expression (33) de- 



rived from the calculations above, we performed a Mon- 
tecarlo simulation, by building numerically the function 
g 2 (R,R,r) for an emitter spectrally diffusing in and out 
of the detection spectral window. We do not take into 
account the subpoissonian statistics of the emitted since 
it can be factorized in the correlation function (see eq. 
|8|. Thus, the simulation gives a direct access to gf n (r) 
and consequently to j3 since 



(35) 



After generating a poissonian stream of photons, we 
assign to each of them a random energy in a gaussian 
distribution. We calculate the probability of the mean 
value of the homogeneous line to jump between the emis- 
sion for the (i — l) th photon and the i th photon. 



7 



The energy of the (i — 1) photon is in the interval: 
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£= 1.7meV 



with SE an infinitesimal energy. 

In the case of a jump of the homogeneous line, the 
energy of the i th photon is in the interval: 
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Figure 6. Calculated half line autocorrelation functions of a 
poissonnian emitter spectrally diffusing, a) for an infinitely 
small homogeneous linewidth, b) for a homogeneous linewidth 
a — 0.3meV. 



Thus, the probability for the homogeneous line to jump 
between the emission of the ith and the (i — l) th photon 
is: 



D jum P (i) = p({^ T dout} | {fJpelin}) = — [l-ezp(-7 c T) 

7c 



(36) 



This probability has already been calculated in the pre- 
vious section (see equation ([29])) . 

7 ou t is the exit rate of the homogeneous line, from I in 
to I ut ilc is the "jump rate" of the homogeneous line. 

For an infinitely small 5E, j c ^ j Q ut and: 



Pjump(i/(i 



!)) = [1 - exp(-j c r) 




If the i photon does not jump, it takes the same 
energy than the (i — l) th photon. 

We then discriminate the photons belonging to ener- 
gies outside the detection area by applying an energy 
condition. They will not be counted in the correlation 
process. 

To calculate the correlation function of the resulting 
photons stream, we compute the delay between the ar- 
rival of each photon and all the other photons of the 
stream. By building the histogram of these delays, one 
obtains the correlation function of the stream [15J. 

The result of the simulation is plotted on fig. [6] a). We 
find again the result explained in the previous section, 
the bunching part of the HLAF takes a value of 2 at zero 
delay. To evaluate j3 from the simulated data, one only 
need to extract g 2 (0). 

To simulate the finite homogeneous linewidth effect, 
we add for each photon an energy shift e. e is a ran- 
dom variable distributed along a lorentzian distribution 
of linewidth a. The effect of this addition is shown on 
fig. |6]b). As expected the bunching is less marked and 
the (3 factor goes from 1 to 0.5 in the case of a finite 
homogeneous linewith of a = 0.3meV and a fluctuation 
amplitude of E = 1.7 meV. 

To compare this simulation with the analytical expres- 
sion, we fix the E parameter and change a. We then re- 
port the values of the f3 factor versus the ratio g in both 
cases (see fig. 0. We can notice that the two methods 
(analytical and Monte Carlo) give the same dependence. 



Figure 7. f3 factor versus the ratio between homogeneous 
linewidth a and fluctuation amplitude S for fixed value of E. 
Red line: analytical model described by eq. [34] Blue dots: 
Monte Carlo simulation result. 



VI. HOMOGENEOUS AND INHOMOGENEOUS 
LINEWIDTH DETERMINATION 

The fit of the experimental data allows to extract, as 
two independent parameters, the correlation time r c and 
/3. As shown in the previous section the latter depends 
on the fluctuation amplitude E and the homogeneous 
linewidth a. To estimate them separately, we can use the 
emission spectrum measured, and by making the same as- 
sumptions as before (the homogeneous line is lorentzian 
and the fluctuation distribution is gaussian) we obtain 
a second equation linking these two parameters since, in 
this case, the emission spectrum has a Voigt profile: 



+oo 

5(<t, E,cj) = J lor(a^).gauss(Ti^ — uj)d^ (37) 



By a double fit, one can find the parameters couple 
(cr, E) which is satisfying these two equations, giving ac- 
cess to the amplitude of the spectral diffusion and the 
homogeneous linewidth separately. 
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VII. CONCLUSION 



The HLAF is not only a function of the correlation 
time of the spectral diffusion. It is dependent on the 
nature of the photon energy statistics. From the calcu- 
lation we performed, one can interpret analytically the 
importance of these different statistics in the observed 
correlation function. By measuring the HLAF and by 
making a temperature dependence of a spectrally diffus- 
ing emitter it is possible to describe how electronic fluc- 
tuations and phonon broadening are evolving, and which 
one is dominating the spectrum. Thus, PCS technique can 
bring all the informations for the full characterization of 
a spectrally diffusing emitter. Indeed, after data treat- 
ment technique presented in the last section of this letter, 
one can obtain separately a (homogeneous linewidth), E 
(fluctuation amplitude) and j c (correlation rate). Other 
existing methods |16j [17] , except the the photon correla- 



tion Fourier spectroscopy method (PCFS) [HI HI] are not 
bringing all these informations at the same time and all 
of them are not adapted in the case of fast fluctuations. 
From an experimentalist point of view, it is important to 
note that the spectral resolution of the PCS experiment 
is limited by the spectrometer, which makes it less accu- 
rate than PCFS , which takes advantage of the Fourier 
transform spectroscopy. However the PCS technique is 
extremely easy to setup and does not require the drastic 
stability conditions demanded by the PCFS technique. 
Since photons correlation experiments require long inte- 
gration time, such sensitivity led to the degradation of 
the time resolution to 20 (is for the PCFS technique, 
which forbids to probe fast fluctuations of the energy 
emitter. In PCS technique, time resolution is only lim- 
ited by photodiodes and can be lowered to around 100 ps. 
This technique is a reliable and accessible technique for 
a complete and fast characterization of a single photon 
emitter in a solid state environment. 
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